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We have carried out thickness dependent first-principles electronic structure calculations on ultra-thin 
films of transition-metal dichalcogenides MX 2 (M=Mo or W; X=S, Se, or Te). When spin-orbit coupling 
(SOC) is included in the computations, monolayer MX 2 thin films display spin-split states around the 
valence band maximum at the Brillouin zone corners with nearly 100% spin polarization. The spins are 
aligned oppositely along out-of-the-plane direction at the K and K' points. For the bilayer films, spin 
polarization of this SOC induced band splitting can be switched on or off by an out-of-the-plane external 
electric field. The spin-polarized states are weakly coupled between the layers in bulk MX 2 compounds with 
small k z dispersion. We confirm a transition from an indirect to direct band gap as the thickness is reduced 
to a monolayer in MoS 2 , in agreement with recent experimental findings. Owing to the presence of a large 
spin- splitting energy and an insulating band gap, MX 2 compounds have great potential for spin/ valley 
electronic applications at room temperature. 

Two-dimensional (2D) materials have been drawing increasing attention because of their potential applica- 
tions in next generation electronics. Graphene 1 , which is currently the most widely investigated 2D material, 
possesses high mobility 2 and displays rich physics 3 " 5 driven by the presence of Dirac cone dispersions near 
the Fermi level at K and K' points at the Brillouin zone corners 6 . However, graphene is a gapless system and the 
spin on the Dirac cone is a pseudospin, limiting its usefulness in semiconductor industry or spintronics device 
technologies. For these reasons, it is clearly important to find new 2D materials, which contain not only a large 
band gap but also support spin-polarized states. 

A number of 2D materials with spin -polarized states have been synthesized. These include, HgTe quantum 
wells and related systems 7,8 , silicene 9 , and transition-metal dichalcogenide (MX 2 ) family 10 " 12 . While HgTe 
quantum well systems and free-standing silicene and related compounds are predicted to be 2D quantum spin 
Hall (QSH) insulators with small band gaps 13 " 15 , MX 2 family harbors much larger band gaps, and a monolayer 
MoS 2 transistor has been reported recently 16 . Radisavljevic et al. used a hafnium oxide gate dielectric to dem- 
onstrate a room-temperature monolayer MoS 2 mobility of at least 200 cm 2 V" 1 s" 1 16 . In monolayer MoS 2 , both 
the valence and conduction band extremes are located at the corners (K-points) of the 2D hexagonal Brillouin 
zone 17 " 19 , resulting in two inequivalent valleys for low energy carriers which are similar to graphene. Bulk MoS 2 
has been demonstrated with a large indirect gap, while its monolayer thin-film exhibits a direct band gap 20,21 . 
Unlike graphene, monolayer MoS 2 lacks spatial inversion symmetry and has a strong spin-orbit coupling (SOC) 
originating from the d orbitals of the heavy transition-metal atoms 19 . It has been shown that due to the presence of 
a large direct band gap, inversion symmetry breaking, and a strong SOC, optical pumping with circularly 
polarized light can achieve a valley polarization of 30% in pristine monolayer MoS 2 22,23 . In order to achieve 
continuous and reversible control of valley- contrasting properties, Wu et al. investigated bilayer MoS 2 by using 
polarization-resolved photoluminescence 24 . In bilayer MoS 2 the circularly polarized photoluminescence can be 
continuously tuned from —15% to 15% as a function of gate voltage, whereas in structurally non-centrosym- 
metric monolayer MoS 2 the photoluminescence polarization is gate independent 24 . Recently, Cheng et al. have 
demonstrated that monolayer MoS 2 is a 2D diluted magnetic semiconductor 25 in which valley polarization can be 
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induced and controlled by magnetic doping 26 . They have also inves- 
tigated the interlayer coupling in ultra thin MoS 2 27 and the Rashba 
effect in MXY (M = Mo, W and X, Y = S, Se, Te) compounds 28 . A 
field-dependent unconventional Hall plateau sequence 29 is predicted 
theoretically in trilayer MoS 2 , and multilayer MoS 2 has been shown 
to be a superconductor at optimal doping 30 . 

Motivated by the intense current interest in 2D transition-metal 
dichalcogenides, the present work undertakes a comprehensive first- 
principles investigation of ultra-thin films of these materials, and 
systematically analyzes their spin-polarization characteristics and 
how these evolve with film thickness, and delineates how different 
chalcogens in Mo and W based films modify their electronic struc- 
tures. When the SOC is included in the computations, band struc- 
tures display spin-split valence bands around the K-points for 
monolayer MX 2 thin-films which do not possess inversion sym- 
metry. The spin of the top valence states is found to be aligned along 
out-of-the-plane direction with the spin pointing up at the K-point 
and down at the K' -point. The inversion symmetry is restored in 
bilayer thin films. Here, we find that when an out-of-the-plane elec- 
tric field E z is applied, the inversion symmetry is removed and the 
spin degeneracy at the K-point is lifted. The spin -splitting energy 
increases with increasing E z and saturates for E z larger than a critical 
field E c . In the high temperature limit, spin polarization for slabs with 
even number of layers is zero, while slabs with odd number of layers 
obey II x decay with increasing film thickness x. Finally, we show that 
the position of the valence band maximum shifts from the K-point to 
the r point as the film thickness increases. In particular, we confirm a 
transition from an indirect to a direct band gap as the thickness is 
reduced to a monolayer in MoS 2 , in agreement with recent experi- 
mental findings 20,21 . 

Computational Details 

MX 2 compounds crystallize in a layered 2H prototype structure with 
space group P6 3 /mmc (Fig. 1), which contains two inverse MX 2 
layers. In each layer, intermediate M atom is sandwiched between 
X atoms, forming strong ionic bonds within a trigonal local structure 
(Fig. 1 (b)). While the intralayer M-X bonding is strong, the inter- 
layer bonding is weak as it arises from van der Waals forces (Fig. 1 
(c)). Unlike the bulk crystal, the space group of monolayer MX 2 is not 
P6 3 /mmc but reduces to P6m2 due to loss of inversion symmetry. 
The inversion symmetry is restored in bilayer thin films: In Fig. 1, we 
used MoS 2 as an example where the inversion center is marked by red 
dot in Fig. 1 (c) lies between two adjacent MoS 2 layers. 

The electronic structures of MX 2 films were computed using the 
projector augmented wave method 31,32 as implemented in the 
VASP 33 " 35 package within the generalized gradient approximation 
(GGA) 36 scheme. The SOC effects were included in a self-consistent 
manner. Slabs with vacuum thickness larger than 15 A were used to 
model the thin films. A15X15X1 Monkhorst-Pack k-point mesh 
was used in the computations. The experimental lattice constants 
were used 10,37,38 . Note that although experimental lattice constants 
are available for all six compounds studied, the experimental atomic 
positions are not available for all these compounds. In order to 
meaningfully unfold trends in electronic structures as a function of 
the number of layers or strain, it is important of course to use a 
consistent theoretical framework in computations. Accordingly, 
relaxed theoretical atomic positions were systematically used in mak- 
ing such comparisons. The atomic positions were relaxed until the 
residual forces were less than 0.005 eV/A. 

Results and Discussion 

Fig. 1 (d) shows the electronic band structure of monolayer MoS 2 . 
Without SOC (red dashed line), a direct band gap of —1.6 eV is 
obtained between the valence band maximum and conduction band 
minimum at the K-point, and all bands are found to be spin degen- 
erate. When the SOC is turned on (black line), the top valence bands 




Figure 1 | (a) Top view of a monolayer MoS 2 thin-film. Gray and yellow 
dots denote Mo and S atom, respectively, (b) Trigonal local structure of 
MoS 2 . (c) Side view of bilayer structure of a MoS 2 thin-film. Red dot marks 
the spatial inversion symmetry point, (d) Band structure of a MoS 2 
monolayer. Black and red dashed lines give results with and without SOC, 
respectively, (e) Spin polarization near the K-point around E F . Color bar 
denotes strength of spin polarization, (f) Spin decay behavior of 
polarization strength near the K-point. 

display a significant splitting (—150 meV) at the K-point due to the 
breaking of spatial inversion symmetry (E(k, f) = E( — k, j)). In sharp 
contrast, the top valence bands are spin degenerate at the T and M 
points. This is because T and M points are time-reversal invariant 
while K-point is not. Focusing on the split band, we find that it 
displays nearly 100% out-of-the-plane spin polarization (Fig. 1 
(e)), which is consistent with previous results 19,39,40 . The two bands 
at K and K' have opposite spin polarizations, obeying the constraint 
of time reversal symmetry (E(k, j) = E( — k, j)). Fig. 1 (f) shows the 
degree of spin polarization for the top valence band near the K- point. 
The bands are nearly 100% spin polarized at the K-point, and decay 
very little as the momentum deviates from the K-point within 0.5 

We emphasize that the presence of a robust 100% spin-polarized 
band at the top of valence bands in an MX 2 monolayer results from 
symmetry constraints of mirror symmetry 40 and the broken inver- 
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sion symmetry of the underlying crystal structure, and therefore, this 
polarization cannot be switched off by an external field. In this con- 
nection, we have investigated MX 2 bilayers in which the inversion 
symmetry is restored with a center of inversion in the structure lying 
between the two layers in the bilayer. Fig. 2(a) shows the band struc- 
ture of a MoS 2 bilayer. Since the inversion symmetry is now restored, 
all bands are spin-degenerate. The energy difference between the top 
of the valence bands at the K and T points, A Kr , becomes smaller. 
Notably, the splitting at the T point is due to the interaction between 
the two MX 2 layers, but the SOC is responsible for the splitting at the 
K-points [Fig. 2(b)]. 

The inversion symmetry of bilayer MX 2 can, however, be removed 
by an out-of-the-plane external field, E z , which lifts the spin degen- 
eracy of states away from the time-reversal-invariant points in the 
Brillouin zone. Since K-points are not time-reversal-invariant points, 
the spin -splitting at the top of the valence bands can be manipulated 
via E z . Fig. 2(b) shows how the bilayer MoS 2 band structure evolves 
with increasing E z . For non-zero E z , there are four spin-polarized 
states in a bilayer, two of which come from the first layer and the 
other two from the second layer in the bilayer system. We define the 
interlayer splitting energy, A inter , to be the energy difference between 
the first-layer spin-up state and the second layer spin-down state. We 
also define the intralayer splitting energy, A intra , to be the energy 
difference between the first- (second-)layer spin-up state and the 
first- (second- )layer spin-down state. At E z = 0.1 eV/A, the states 
from the second layer are slightly lower than those from the first 
layer in each split branch, making A inter less than A intra , in agreement 
with previous results 41 " 43 . [The two topmost valence states are from 
different layers.] When E z is larger than a critical value E c , A inter 
becomes larger than A intra , making the two topmost valence states 
to come from the same layer. Fig. 2 (c) shows A inter and A intra as a 
function of E z . We see that A inter increases linearly as E z increases, but 
Aintra remains nearly constant. A inter is due to the potential difference 
between two layers caused by E z . On the other hand, A intra is the spin- 




Figure 2 | (a) Band structure of bilayer MoS 2 . Black and red dashed lines 
give results at E z = 0 and E z = 0.3 (eV/A), respectively, (b) Evolution of 
band structure near the K-point (red dashed box in Fig. 2(a)) for varying 
external electric field, E z . Color scale gives the degree of spin polarization. 
Contributions from the first and second layers in MoS 2 bilayer are marked 
with numbers 1 and 2. The unit of electric field is eV/A. (c) Spin- splitting 
energy at the K-point as a function of E z . Blue line gives the intralayer 
splitting energy, A intra , While the red line gives the interlayer splitting, 

Ainter- 



splitting resulting from the effect of SOC within each layer, and it 
therefore depends weakly on E z . There is a crossover between A inter 
and 

Aintra for 0.1 < E z <C 0.2. This indicates that the maximum of the 
spin-splitting energy is about 150 meV (solid line in Fig 2(c)), which 
is comparable to the intrinsic value of the monolayer spin -splitting 
(Aintra)- A inter under high E z (E z > E c ), red dashed line in Fig 2(c), is 
not relevant to the spin- splitting at the two topmost valence bands, 
both of which now primarily come from the same layer. Note that the 
main effect of E z is to induce a significant band- splitting with little 
change in the degree of spin polarization. The spin direction reverses 
when the direction of the E z reverses, providing an effective pathway 
for manipulating spin states. These spin-polarized states could be 
observed via spin sensitive angle -resolved photoemission experi- 
ments, important matrix element effects in ARPES and other highly 
resolved spectroscopies 44 notwithstanding, and should be accessible 
to transport experiments in hole-doped samples. 

Figure 3(b) shows the band structure of bulk MoS 2 at selected k z 
values. The two well-separated branches can still be seen at the K 
point because the SOC induced splitting is much larger than the k z 
dispersion resulting from interlayer coupling. The SOC splitting is 
about 150 meV, which is comparable to the value of A intra in the 
bilayer or the spin -splitting energy in the monolayer, while the inter- 
layer coupling at the K-point is estimated to be only about 25 meV. 
Note that when the number of layers is finite, the number of states in 
each branch is the same as the number of MX 2 layers. 

Concerning the spin-polarization for the top branch of the valence 
bands at the K-point, it is interesting to consider the degree of spin 
polarization in the high temperature limit when states in the top 
branch will become accessible due to thermal fluctuations. Since a 
multilayer with an even number of layers is inversion symmetric, 
there will be no net spin polarization in a branch containing an even 
number of states. Adding one MX 2 layer to create a multilayer with 
an odd number of layers, however, adds a fully spin polarized state 




Number of layers 



(b) kz=0.0 (2n/c) — -kz=0.1(27i/c) 




Figure 3 | (a) Average spin polarization of the top branch of valence bands 
as a function of the number of MX 2 layers, (b) Band structure of bulk MoS 2 
at selected Rvalues (in unit of 27i/c). Solid and dashed lines of various 
colors refer to different Rvalues (see legend). 
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into the branch. Despite the mixing of other states through the weak 
interlayer coupling, the average out-of-the-plane spin polarization 
for the whole branch becomes nonzero. The degree of spin polariza- 
tion for the top branch is inversely proportional to the number of 
layer, x. As shown in Fig. 3(a), the average spin polarization, < S z >, 
of the top branch at the K-point as a function of thickness in the high 
temperature limit can be expressed as 



0 x = 2,4,6... 
fx x= 1,3,5... 



Note that the states at K' possess opposite spin polarization. 
Therefore, f(x) will take minus sign for the top branch at K'. 

We comment further on the nature of valence bands at the T point. 
Unlike the K-point, since T is a time-reversal-invariant point, states 
at r are spin degenerate for any number of layers. With increasing 
number of layers, as more states interact through interlayer coupling, 
the bandwidth increases, and the top of the valence band at T moves 
to higher energy. In MoS 2 , the interlayer coupling at T is —0.25 eV, 
an order of magnitude larger than that for top valence states at the K- 
point. This drastic difference in the strength of interlayer coupling is 
due to the orbital character of the relevant states: The top valence 
states at T are primarily associated with out-of-the-plane d 3z 2_ r 2 
orbital, while states at the K-points arise from in-plane d xy and 
d x 2_ y i orbitals. 

Accessibility of spin polarized valence states for spintronics appli- 
cations will be controlled by the size of SOC spin-splitting energy 
A intra and the energy difference between the top valence states at K 
and r, A KY = E K - E r . Accordingly, Fig. 4(a) compares these key 
quantities in various MX 2 materials. A intra is seen to increase with 
increasing atomic number as expected since the strength of SOC is 
larger in heavier atoms. A KY , on the other hand, does not show any 
simple trend. Interestingly, the value of A KY correlates with that of X- 
M-Xbond angle (solid green line). This indicates that the crystal field 
and M-X bonding, which determine the relative positions of various 
orbitals, also play an important role in determining A Kr . 

In order to help guide search for suitable substrates, we have 
explored the effect of strain on A intra and A KY . Fig. 4(b) shows A KY 
as function of strain for monolayer, bilayer, and bulk MoS 2 . A KY is 
seen to increase (decrease) as the lattice constant decreases 
(increases). This trend is consistent with our observation above that 
A KY correlates with the X-M-X bond angle in that thin films with 
larger in-plane lattice constant possess smaller X-M-X bond angles. 
Interestingly, A KY is very small for bilayer MoS 2 , and only a 1% lattice 
constant increase changes its sign, adding to the possibilities for 
control of spin-polarized states via an external electric field. 
Notably, the small value of A^r in bilayer can change sign if experi- 
mental atomic positions are used in the computations instead of self- 
consistently optimized theoretical positions 41 " 43 . This discrepancy 
between experiment and theory is due to overestimation of the 
bond-length and bond- angle in GGA, which is a typical error of 
the GGA exchange-correlation functional 45 . Despite these uncertain- 
ties in absolute band positions, the trends shown are fairly robust to 
whether theoretical or experimental atomic positions are used con- 
sistently in the calculations. We suggest that substrates should be 
chosen with smaller lattice constants so as to have larger A KY values. 
Note also that unlike A KY , which is sensitive to in-plane strain and 
slab thickness, A intra is found to be almost constant (not shown) 
when strain, E z , or thickness is varied. 

Concerning the nature of the band gap and whether it is direct or 
indirect, we note that all MX 2 considered here have valence band 
maximum at T for bulk and at K for monolayer. Take MoS 2 as an 
example. Fig. 4(c) shows that A KY decreases as slab thickness 
increases. As shown in Fig. 3(b), the valence band maximum is at 
r and conduction band minimum is away from T along the T-K 
direction, so that the bulk band gap is indirect. For monolayer MoS 2 
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shown in Fig. 1 (d), on the other hand, both the conduction minimum 
and valence bands maximum shift to the K-point, making a direct 
band gap. Our results are consistent with the fact that bulk MoS 2 has 
been demonstrated to support a large indirect-gap, which crosses 
over from indirect to a direct gap for a monolayer 20 ' 21 . Such an indir- 
ect to direct band gap transition as the thickness is reduced has also 
been observed directly via ARPES spectroscopy for MoSe 2 46 . 

Conclusions 

We have systematically investigated electronic structures of mono- 
layer and bilayer films of transition-metal dichalcogenides MX 2 
(M=Mo or W, X=S, Se, or Te) via first-principles calculations, 
and discuss how these electronic structures evolve in the bulk limit. 
In the monolayer case, the spin -splitting due to SOC occurs at the 
valence band maximum, resulting in two non-degenerate nearly 
100% spin-polarized states with opposite out-of-the-plane polariza- 
tions at the K and K' symmetry points in the Brillouin zone. We show 
that the SOC induced spin splitting can be tuned by an out-of-the- 
plane external electric field for bilayer MX 2 thin films. Our theor- 
etical results concerning the nature of the band gap, whether it is 
direct or indirect, and its evolution with layer thickness are in accord 
with recent experimental findings in several MX 2 materials. Our 
study provides insight into the degree of spin polarization of the 
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valence bands in ultra-thin transition-metal dichalcogenide films 
and their viability for spintronics applications. 
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